CHOQUET ORDER FOR SPECTRA OF HIGHER LAME 
OPERATORS AND ORTHOGONAL POLYNOMIALS 



JULIUS BORCEA 

Abstract. We establish a hierarchy of weighted majorization relations for the 
singularities of generalized Lame equations and the zeros of their Van Vleck and 
Hcine-Stieltjes polynomials as well as for multiparameter spectral polynomials 
of higher Lame operators. These relations translate into natural dilation and 
subordination properties in the Choquet order for certain probability measures 
associated with the aforementioned polynomials. As a consequence we obtain 
new inequalities for the moments and logarithmic potentials of the correspond- 
ing root-counting measures and their weak-* limits in the semi-classical and 
various thermodynamic asymptotic regimes. We also prove analogous results 
for systems of orthogonal polynomials such as Jacobi polynomials. 



1. Introduction 

The generalized Lame equation in algebraic form is the second order differential 
equation 

Q 2 (z)y"(z) + Q 1 (z)y'(z) + Q () (z)y(z) = 0, (1.1) 

where Q 2 , Qi,Qo G C[z] with degQ 2 = P, degQi = p — 1, degC? < P - 2. Par- 
ticularly important cases are p = 2 and p = 3, which correspond to the hyper- 
geometric differential equation and Heun's equation, respectively (cf. [23 ). The 
classical Heine-Stieltjes multiparameter spectral problem deals with so-called Lame 
solutions of the first kind (of given degree and type) to equation (jl.ip and may be 
formulated as follows: given Q2 (z) , Qi (z) as above and n G N find a polynomial 
V{z) of degree at most p — 2 and a polynomial S{z) of degree n such that (jl.ip holds 
for Qo(z) = V(z) and y(z) = S(z). If such V(z) and S(z) exist we say that (jl.ip is 
n-solvable. A generalized Lame equation is solvable if it is n-solvable for all n G N. 
The coefficients V(z) are called Van Vleck polynomials and the corresponding so- 
lutions S(z) are known as Heine-Stieltjes polynomials. These two classes are also 
referred to as Lame polynomials or generalized spectral polynomials for (jl.ip . 

There are several known sufficient conditions for the solvability of equation (|l.ip . 
For instance, Heine 18] proved that for any n G N there exist at most 

different Van Vleck polynomials V(z) for which (jl.ip has a polynomial solution 
y(z) = S(z) of degree n. Heine's text is written in a traditional XlXth century style 
German and the exact statements it contains seem to have created some confusion 
(cf. [H]). Szego [551 §6-8] quotes this result and adds that "Heine asserts that, in 
general, there are exactly cr(n) determinations of this kind". As explained in [71 ITU]. 
Heine actually showed that if the coefficients of Qi(z) and Qi{z) are algebraically 
independent - that is, these coefficients satisfy no algebraic equation with integer 
coefficients - then (jl.ip is solvable. Moreover, if this is the case then for any neN 
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there exist exactly a(n) different Van Vleck polynomials V(z) of degree p — 2 and 
the same number of corresponding monic Heine-Stieltjes polynomials S(z) of degree 
n. An explicit characterization of the exceptional cases when this number is strictly 
less than a(n) seems to be lacking for the moment [22 . In the general case, Heine's 
arguments imply that (jl.ip is n-solvable for all sufficiently large n [TP] , 

The solvability of (jl.ip has been established under various other assumptions, 
most notably when Q2(z) and Q\{z) have strictly interlacing real zeros and the 
leading coefficient of Qi(z) is positive. This case is particularly interesting from a 
physical point of view and has attracted a lot of attention in recent years. Indeed, 
differential equations of the form (jl.ip whose coefficients satisfy the above condition 
arise naturally when separating variables in the Laplace equation in spherical co- 
ordinates and yield important examples of quantum completely integrable systems 
such as generalized (real or complex) Gaudin spin chains [UJ [T51 [J_3] • A fundamen- 
tal result of Stieltjes [25] - also known as the Heine-Stieltjes theorem [26] - asserts 
that if Q2(z) and Q\{z) have strictly interlacing real zeros and Qi(z) has positive 
leading coefficient then for each n £ N there are exactly a(n) different Van Vleck 
polynomials V(z) of degree p — 2 and the same number of corresponding monic 
Heine-Stieltjes polynomials S(z) of degree n. The latter are given by all possible 
ways of distributing the zeros of S(z) in the p — 1 open intervals defined by the 
zeros of Q2(z). Stieltjes actually showed that the zeros of S(z) are the coordinates 
of the equilibrium points of a certain electrostatic potential. Similar results have 
recently been obtained in cases when Q2{z) has all real zeros and the residues in 
the partial fractional decomposition of Qi(z)Q2(z) _1 have mixed signs [TBI, I17j . 

Let us assume that Q2{z) and Qi(z) are such that 



Note that if Q2{z) and Q\(z) are as above and equation (jl.ip is solvable then 
any Van Vleck polynomial is of degree exactly p — 2. Polya [51] - and Klein and 
Bocher before him (cf. [TBJ) - showed that in this case the zeros of all Van Vleck 
and Heine-Stieltjes polynomials lie in the convex hull of £i, . . . , £ p . Extensions of 
this Gauss-Lucas type theorem to cases when the residues a.; are not necessarily 
positive real numbers as well as various other results on the location of zeros of 
Lame polynomials have since been obtained [JJ [57] ■ In this paper we show 
that much more is actually true. Namely, if (jl.2p holds then the zeros of any Van 
Vleck polynomial together with those of a corresponding Heine-Stieltjes polynomial 
and the zeros of Q2(z) satisfy a weighted majorization relation in the sense of [3] 
(see §2). This amounts to a dilation property - equivalently, a subordination re- 
lation in the Choquet order - for certain probability measures associated with the 
generalized spectral polynomials and the singularities of equation (jl.ip . A precise 
statement of this result is given in Theorem [3J below. As a consequence we obtain 
new inequalities for the moments and logarithmic potentials associated with the 
root-counting measures of Lame polynomials and we establish similar properties 
in the thermodynamic (p — > oo) and semi-classical (n — > oo) asymptotic regimes 
(Corollaries [JH3J). These results hold in the greatest possible generality and require 
no additional assumptions besides (|1.2p . Therefore, Theorem [3] and Corollary [JJ 
apply whenever equation (jl.ip is n-solvable while Corollaries [5H3J make sense in all 
cases when (jl.ip is solvable and the considered limits exist (see $3j for several con- 
crete examples) . In the special case when Q £ R, 1 < i < p, our results are a natural 
complement to those of [7J [JJJ [J5J [J3] [53] dealing with asymptotic distributions, 
limiting level-spacings and mean densities of zeros of Lame polynomials. 




(1.2) 



where Q G C and a\ > 0, 1 < I < p. 
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Various extensions of the Hcinc-Stieltjes multiparameter spectral problem to 
higher order linear ordinary differential operators with polynomial coefficients have 
been studied in 0110]. In particular, if Qi(z) and Q\{z) are as in (II. 2j) and k>2 
then one may consider an operator of the form 

5(z) = ° 2(z) £ + Ql(z) £^T- (!-3) 

As in jlOj , we call a higher order generalized Lame operator or a higher Lame 
operator for short, provided that its Fuchs index r := p — k is non-negative. If 
r = then £>(z) is a so-called hypergeometric type operator. Such operators and 
their polynomial eigenfunctions have important applications to the study of the 
Bochner-Krall problem and exactly solvable models (see, e.g., [3 HUH] and refer- 
ences therein). The multiparameter spectral problem for a higher Lame operator 
t)(z) is as follows: given n 6 N find a polynomial V(z) of degree at most r such 
that the equation 

Hz)y(z) + V(z)y(z) = (1.4) 

has a polynomial solution y(z) = S(z) of degree n. One can then define the notions 
of n- solvability, solvability, higher Van Vleck and Heine- Stieltjes polynomials - that 
is, higher spectral polynomials or Lame polynomials - corresponding to (|I.4p by 
analogy with the terminology used for Several sufficient conditions for the 

solvability of (|I.4p that extend those of Heine for (jl.ip were recently obtained in [10] 
(see §3). We show that whenever equation (|1.4[) is solvable its singularities and the 
zeros of all corresponding higher Lame polynomials satisfy weighted majorization 
relations (Theorem 0]) and we establish natural analogs of Corollaries HH3] for the 
higher order case (Corollaries [HE]). 

Our methods also yield interesting applications of the Choquet order / weighted 
majorization to the theory of orthogonal polynomials. In particular, we prove ap- 
propriate versions of the aforementioned results for classical orthogonal polynomials 
such as (ultraspherical) Gegenbauer polynomials, (associated) Legendre polynomi- 
als, Chebyshev polynomials and indeed any family of Jacobi polynomials (see p.3p . 

This paper is organized as follows. In $2]we recall the notion of weighted multi- 
variate majorization from [5] as well as the definition and properties of the Choquet 
order for non- negative Radon measures. We state and prove our main results in fJ3J 
In 2] we give several generalizations and discuss some related problems. 



2. Weighted Majorization and the Choquet Order 

The majorization preorder on n-tuples of real numbers - also known as the strong 
spectral order, vector majorization or classical majorization - essentially quantifies 
the intuitive notion that the components of a real n-vector are less spread out 
than the components of another such vector. Several matrix versions of this notion 
have been proposed and studied in various contexts |21j . A weighted multivariate 
extension of both classical and matrix majorization was introduced in [3j. In the 
special case of complex n-vectors the definition of loc. cit. is as follows. For m G N 
set 

Am = < a = (ai, . . . , a m ) 




X m = C m xA m , X=|JX m . 

Denote by conv(Jl) the (closed) convex hull of a (bounded) set SI C C and by 
X T the transpose of a (row) vector X ~ (xi, . . . , x m ) £ C m . We frequently write 
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conv(X) for conv({a;i, . . . , x m }). Let n be the set of all row stochastic m x n 
matrices. 

Definition 1. The pair (X, a) £ X m is said to be weightily majorized by the pair 
(Y, b) £ X„ , denoted (X, a) -<! ( Y, b) , if there exists a matrix R £ MJ^ n such that 

X T = RY T and b = aft, 

where X T and Y T are obtained by some (and then any) ordering of the coordinates 
of X T and Y T , respectively. 

Remark 1. Note that if (X, a) -< (Y, b) then X £ conv(Y) m and the a-barycenter 
of X must coincide with the b-barycenter of Y, that is, J2iLi a i x i — ^jVj- 
Moreover, it is clear that the weighted majorization relation is both reflexive and 
transitive, which makes it a preorder on X. One can also show that for every m £ N 
this preorder induces a partial order on the orbit space C m /E m , where E m is the 
symmetric group on m elements. 

The following characterization of the weighted majorization relation may be 
found in [3l Theorem 1]. 

Theorem 1. Let (X, a) £ X m and (Y,b) £ X„, where X = (x\, . . . , x m ) £ C m . 
a = (<xi, . . . ,a m ) € An, Y = (yi, . . . ,y n ) € C" and b = . . . ,b n ) € Ai- ^fte 
following conditions are equivalent: 

(i) /or any (continuous) convex function f : C — > R one ftas 

m n 

^ai/(xi) < ^bjfiyj); 

«=1 3=1 

(ii) i/ie relation (X, a) -< (y, b) holds. 

Remark 2. If (X, a) -< (Y, b) then the inequality in Theorem [1] (i) holds for every 
convex function / defined on conv(Y). 

There is a natural connection between the weighted multivariate majorization 
relation and the Choquet order for non-negative Radon measures. The latter has 
been studied in the general context of locally convex separable topological vector 
spaces in e.g. [H] and subsequent papers. For measures defined on compact subsets 
of the complex plane the Choquet order and the main results of op. cit. may be 
described as follows. Let K be a convex compact subset of C, denote by C(K) the 
space of real continuous functions on K and let V(K) the subset of C(K) consisting 
of convex functions. If /z is a non-negative Radon measure on K and / is a function 
on K one defines fi(f) = J K f(y)dfi(y). The mass of \i is therefore — f„ dfi(y) 
and if > then the barycenter of \i is the point r(/i) := /i(I) _1 J„ yd^,{y). 

Definition 2. Given two non-negative Radon measures fj, and v on K one says 
that v dominates \x in the Choquet order or that v is a dilation of /i, denoted /i -< is, 
ifM/)<K/)forany/eP(X). 

The use of the term "dilation" in Definition [5] is motivated by Definition [3] and 
Theorem [2] (iii) below. To formulate this result we need a few more concepts and 
notations. Let M.{K) be the set of all probability measures with supp/i C K . Note 
that if /i £ A4(K) then its barycenter r(/i) lies in K . 

Remark 3. As is well known, the set Ai(K) equipped with the weak-* topology is 
a sequentially compact Hausdorff space. This will allow us to choose a convergent 
subsequence from any sequence of measures belonging to A4(K). 
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Definition 3. A dilation on A is a weakly Borel measurable M. ( A)-valued function 
on A that inverts the barycenter mapping. In other words, a map T : K — > A4(K), 
x i— ► T x , is a dilation on A if r{T x ) — x for all x G A and the real- valued function 
on A given byii-> T x (f) is borelian for any / G C(K). 

If T is a dilation on if then for any non-negative Radon measure fi on A one 
can define a new such measure v :— T(p) by setting 

K/) - / T x (f)dn(x), f e C(A). (2.2) 

It is not difficult to show that the real- valued function on A given by x i— > T x {f) is 
borelian and bounded whenever / is a bounded borelian real-valued function on A' 
and that (|2.2|) actually holds for all such functions (cf. [2]). The main results of 
loc. cit. provide various descriptions of the Choquet order in a quite general setting. 
In the case discussed above these may be summarized as follows (see also [16]1- 

Theorem 2. If A is a convex compact subset of C and fj,,u € M{K) then the 
following conditions are equivalent: 

(i) \i -< v; 

(ii) for every convex combination \x = X)"=i ^ifM (li G A4(K), 1 < i < n, 
there exists a corresponding convex combination v = ^2™ = i^iVi such that 
Vi G M{K) and r{v{) = r(/ii), 1 < i < n; 

(iii) v — T(/i) 7 where T is a dilation on A. 

Remark 4. If the conditions in Theorem[2]hold then supp(/i) C conv(supp(i/)). This 
may be viewed as a "continuous" version of the corresponding result for weighted 
majorization (cf. Remark [1]). 

3. Main Results and Proofs 

Given a complex polynomial P of degree d > 1 we let Z(P) be the d-tuple (or 
multiset) consisting of the zeros of P, where it is understood that each zero occurs 
as many times as its multiplicity. In particular, |i?(P)| — d. To P we associate its 
root-counting measure, namely the (finite) real probability measure given by 

v P = \z{p)\-^ s <> 

Cez(P) 

where 8q is the Dirac measure supported at £. The symbol V is used below for the 
concatenation operation, that is, if [x\, . . . , x m ) G C m and (yi, . . . , y n ) G C n then 

(xi, . . . , x m ) V {yi, y n ) = (x x , . . . , x m , yi, ■ ■ ■ , y n ) G C m+ ™, 

and the "all ones" vector is denoted by l m = (1, . . . , 1) G K m . 

3.1. Generalized Lame operators. Suppose that n > 2 is an integer such 
that (|1.1|) is n-solvable and that Q2(z), Q\{z) satisfy (|1.2p . Let S(z) be a Heine- 
Stieltjes polynomial of degree n corresponding to a Van Vleck polynomial V(z), so 
that degV = p — 2 (cf. EJTJ) . Let a = a(n,p) := n — 1 + ^? =1 a; and define the 
following weight vectors: 

a— a 1 b— a + n—1 

(p — l)a + n ~ 1 v 2l n\(p — l)a + n — 1] 71 ' 

/ a — ai a — a p 

\(p— l)a + n — 1 ' ' (p — l)a + n — 1 
Recall (|1.2[) and note that a > 1, c G while a V b G A n + P -2- Finally, set 

Z(V) = (v 1 ,...,v p - 2 ), Z(S) = ( Sl ,...,s n ). (3.2) 
We can now state our first main result. 
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Theorem 3. With the above notations and assumptions the inequality 
p-2 r / i\/ 1 M " p 

£/(»,) + i- (i- i) (i- ±) £/(.,) ^ '-3 /«) <«) 

fto/ds /or anj/ convex function f : C — > M and i/ equality occurs in (|3.3|l /or some 
strictly convex function f then the zeros of Q 2 must be collinear. Equivalently, 

(Z(y)VZ(S),aVb) -< (Z(Qa),c). 

T/ims i/iere exists a matrix R = R(n,p) G M,^ s +p _ 2 p swc/i iftai 

(Z{V) V Z(S*)) T = RZ(Q 2 ) T and c = (a V b)i?. 

Remark 5. As pointed out in EJTJ the only requirements for Theorem [3] are that 
p.2p holds and equation is n-solvable. For instance, Stieltjes' theorem shows 
that (jl.ip is always solvable if Z(Q 2 ) C M while Heine's result [TH] asserts that the 
same is true whenever Q 2 (z) and Q\{z) are algebraically independent. 

Note that in particular Theorem [3] immediately implies the Polya-Klein-Bocher 
result mentioned in SjlJ namely Z(V) U Z(S) C conv(Z(Q2)) (cf. Remark[T]). Now 
given a compact set if C C and fi € A^(if) let 



|w| m d/i(w), r7i ez + , 

denote the moments of p. As is well known, the logarithmic potential of p 

U»(z) = J \og\z - w\dp[w) 

is subharmonic in C and U fJ '(z) — — oo for every atom z of p. 

Clearly, (|3.3p may be reformulated in terms of the Choquet order for atomic 
probability measures with finite point spectrum: 

Corollary 1. In the situation of the preceding theorem one has 
(p — 2)a a + 7i — 1 

{p - l)a + n - + (p-l)a + 77-l Ms * ^ ' 

where fi v and p s are the root- counting measures of V and S , respectively, while 
p Q £ A4(conv(Z(Q 2 ))) is defined by 

su PP (/i Q2 ) = Z{Q 2 ) = {0}f =1 and p Q2 ({0}) = ( p _ ") a _ 1 

for 1 < I < p. In particular, 

(P~2)a (m) a + n-l (m) < (m) 

(p- l)a + n- r v (p- l)a + n- I s ° 2 
/or aZZ m £ Z + and 

(P - 2)a U-v (z) + , U X n ~ l , U"s (z) > U*"* (z) 



(p — l)a + n — 1 (p — l)a + 77 — 1 

whenever z € C \ conv(Z(Q2))- 

In the semi-classical asymptotic regime (77—5-00) Theorem [3] yields: 

Corollary 2. Assume that (|1.2[) ftoWs and i/iai is solvable. Let {S , n (-z)}neN 
a sequence of monic Heine- Stieltjes polynomials such that deg i5„ = n, n € N, and 
{V n (z)}neN fre a corresponding sequence of Van Vleck polynomials with degV^ 
fi 1 normalized so that each V n is monic. Then 

P-2 2 
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where /i* = lim^ 9lwoo \i and /i* = hm^gn^oo ^ s „ / or 071 appropriately chosen 
AcN. Equivalently, there exists a dilation T on conv(Z(Q 2 )) s^c/i ttffli PMq 2 — 
T((p — 2)/i* + 2/i* ) . In particular, 

P~ 2 *< m > , 2 *< m ' ^ -Cm) r- n? 



and 



p " p J Q2 



P 2 -U^v(z) + -U^s(z)>U^Q2(z) 



p p 
for any z £ C \ conv(Z((2 2 )). 

Finally, we may also let p — > oo and consider various so-called thermodynamic 
asymptotic regimes (cf., e.g., [ITJ EH] ) ■ In this case we get the following result. 

Corollary 3. Let {Q2.p(z)} p x L 2 an d {0i,p( z )}p^=2 be ^ wo sequences of polynomials 
such that for all p > 2 the pair (Q2, p (z),Qi tP (z)) satisfies (|1.2[) and the corre- 
sponding equation is solvable. Assume further that K C C is a compact set 
with Z(Q2,p) C K, p > 2, and that {S p>n (z)} n ^, respectively {V p , n (z)} ne ?i, is a 
sequence of monic Heine- Stieltjes polynomials, respectively normalized Van Vleck 
polynomials, associated with the resulting system of equations such that deg S PjTl = 
n, deg Vp, n = p — 2 and V p>n is monic for all p > 2 and n € N. Then 

*l L v *1 1 q 2 ' 

where = lim^^ ^ Vp n(p) and „^ Q2 = lim^^ fj,^ for an appropriately 
chosen T C N. Hence there exists a dilation T on K such that *p. Q — T^fx v j. In 
particular, 

*^<^™\ mez + , 

while 

U^v{ z ) > [T m q 2 ( z ) 

whenever z£C \ K . 

3.2. Higher Lame operators. Let now k > 2 be a fixed integer and consider an 
order k generalized Lame operator D(z) with Fuchs index r :— p— k as in (| 1 - 3[) and 
the corresponding multiparameter spectral problem (|1.4p . Assume that the latter 
is n-solvable and that (S(z),V(z)) is a pair of (higher) spectral polynomials with 
deg S = n and degV = r satisfying ()1.4|) . Let a k = at(n,p, k) := n — k + l + Y^—i a i 
and define the following weight vectors: 



_ «fc _ (fc - l)q fc + n - k + 1 

(p — + n. — fc + 1 r ' n[(p — l)o!fc + n — k + 1] 



Cfe 



a fe - ai a k 



(3.4) 



(p — + n — k + 1 (p — l)ci!fc + n — fc + 1 ( 

One clearly has > 1, € „4 P , a^ V b^ e Ai+r, a 2 = a, b 2 = b and c 2 = c. 
Since in this case deg V = r we adapt notation (|3.2p to the current situation simply 
by setting Z(V) = (vi, . . . ,v r ). 

The analog of Theorem [3] for higher Lame operators reads as follows. 



Theorem 4. Under the above assumptions the inequality 



2—1 L X 7 X ' J J — l l — l X ' 

ZioZds for any convex function f : C — > K and i/ equality occurs in (|3.5|) /or some 
strictly convex function f then the Ci 's must be collinear. Equivalently, 

(Z(V) V Z(5), a fe V b fc ) -< (Z(Q 2 ), c fe ) . 
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Thus there exists a matrix R k = R(n,p, k) G NS.™ +rp such that 

(Z(V) V Z(S)) T = R k Z(Q 2 ) T and c k = (a* V b k )R k . 

Remark 6. The proof of Theorem[4] actually yields an inequality stronger than ()3.5p 
involving all four zero sets Z(V), Z(S), Z(S*( fe - 1 )) and Z(Q 2 ) (see (ETTUll below). 

Remark 7. Theorem |4] applies to all situations when Q 2 (z), Qi(z) satisfy (|1.2p and 
equation (|1.4p is n-solvable. By [TUl Theorem 5] this is always true for all sufficiently 
large n. Moreover, it was shown in op. cit. that (| 1 .4|) is 7i-solvable for any n € N in 
each of the following cases: (i) Q 2 (z) and Q\{z) are algebraically independent (ii) 
2{Q 2 ) C K (hi) D(z) is a hyperbolicity preserving operator (HPO for short), i.e., it 
maps polynomials with all real zeros to polynomials with all real zeros. A complete 
classification of all HPOs was recently obtained in [5] (see also [5]). 

Natural extensions of Corollaries [T] EH to higher Lame operators are as follows. 

Corollary 4. In the situation of Theorem [7] one has 

(p - k)a k (k - l)a k + n — 1 



+ ( n _ 1 \ n , _L r, _ 1 ^ S ^ ' 



(p - l)a k v {p - l)a k + n — 1 

where (i v and /i s are i/ie root-counting measures of V and S , respectively, while 
jl Q € A^(conv(Z(Q2))) is defined by 

supp(/i Q J = Z(Q 2 ) = {OILi and Mq 2 ({0}) = 



(p - l)afe 



/or 1 < i < p. 7n particular, 



(p - k)a k (m) (fc - l)Qfe + n - 1 (m) < (m) 
(p- l)a fc +n - V v {p-l)a k +n-V s ~^ Q ^ 



for all m € Z + and 



ip-k)a k (fc-l)a fc + n -!,.„, ,,. 



-ET^v (z) + i -fJLJ -C^s ( z ) > (z) 



(p - l)a fe + n - 1 (p - l)ajfe + n - 1 

whenever z 6 C \ conv(Z(Q2))- 

Corollary 5. Assume that (|1.2[) ZioWs and £/ia£ (|1.4p is solvable. Let {S n (z)} ni zjq be 
a sequence ofmonic higher Heine- Stieltjes polynomials such that degS",! = n, n € N, 
and {Ki(-z)}n6N be a corresponding sequence of higher Van Vleck polynomials with 
deg V n = p — 2 normalized so that each V n is monic. Then 

P-k » fc * 

M v + "Ms ^ Mq 2 , 

p P 

where /i* = lim^g^^^ /i^ and /i* = lim^ 3n ^ 00 /or an appropriately chosen 
AcN. Equivalently, there exists a dilation T on conv(Z(Q 2 )) suc/i iftai PM Q2 = 
T((p — fc)p* + fc^o) • ^ n particular, 

Mv +-M S </4,» meZ+, 

p p Q2 

and 

-U^v (z) + -U^s ( z ) > ( z ) 
P P 

for any z 6 C \ conv(Z(Q2))- 
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Corollary 6. Let {Q2, P (z)} r ^ =2 an d {Qi, P i z )} p x L2 be two sequences of polynomials 
such that for all p> 2 the pair (Q2, P (z), Qi, P (z)) satisfies (|1.2p and the correspond- 
ing higher Lame equation (|1.4p is solvable. Assume further that K C C is a compact 
set with Z(Qi >p ) C K, p> 2, and that {<Sp,n(.z)}neN> respectively {V^ jn (z)} n6 N, is a 
sequence of monic higher Heine- Stieltjes polynomials, respectively normalized higher 
Van Vleck polynomials, associated with the resulting system of higher Lame equa- 
tions such that deg S p ^ n = n, deg V P , n = p — 2 and V PtU is monic for all p > 2 and 
neE Then 

*Mv ~^ *1 1 q 2 ' 

where *^ v . = lim^^ H Vp n(p) and = lim^^ fi Q2 p for an appropriately 

chosen fcN. Hence there exists a dilation T on K such that *M Q2 — ^ n 



particular, 
while 

whenever z G C \ K . 



U' fi v(z) > [T M Q 2 ( z ) 



3.3. Orthogonal polynomials. As we shall now explain, the above results have 
interesting yet apparently unknown analogs for important classes of orthogonal 
polynomials such as Jacobi polynomials. Recall that the latter are defined by 

pt \z) = ^-(i - Z )- a (i + z)-p— [(i - z)" + «(i + zr+p] , 

where a > — 1, [3 > — 1. For special values of the parameters a and (3 one gets 
(up to a normalizing factor) all the other classical Jacobi-like polynomials includ- 
ing (ultraspherical) Gegenbauer polynomials, (associated) Legendre polynomials 
and Chebyshev polynomials of the first or second kind. The relative location and 
asymptotic behaviour of the zeros of Jacobi polynomials have been of permanent 
interest in view of their important role as nodes of Gaussian quadrature formulae 
and their nice electrostatic interpretation (Bethe Ansatz) [26] . 
We prove the following result. 

Theorem 5. Let Z(P^ a ' P) ) = {Cn,i}?=i be the zero set of P^' P) . Then 

IX-f,, W (n + f3)f(l) + (n + a)f(-l) 
»^ AWS 2n + a + (3 

i—l 

for any convex function f : [—1, 1] — > M. 

Remark 8. The proof of Theorem [5] yields in fact an even stronger relation that 
involves both zero sets Z(P^ P) ) and z(p j 1 q ' /3) ') , see ([3~TTj) in gUbelow. 

Remark 9. Note that if /i {a ^ denotes the root-counting measure of Pt'^ then 
Theorem [5] may be restated in terms of the Choquet order simply as 

n + (3 n + a 

^i"-* 3 * 2n + a + P 1 + 2n + a + /3 
In particular, by letting n — > oo we get 

+ (3.6) 

where /x* ^ := lim*^^ [i (q 0) is the *-limiting distribution of the zeros of P^'^\ 
As is well known (see, e.g., [3S]) the latter is the (uniform) arcsine distribution and 
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thus (|3.6[) may be rewritten as 



i r m __ dz K /(D + /(-i) 



for any convex function / : [— 1, 1] — > R. Although elementary, the above inequality 
is not completely obvious; arguably the most direct way of proving it is to note 
that it actually amounts to showing that 

i r 1 \z-c\ 
W-i VT 

for any c £ R, which is a trivial exercise. 



f 1 \z — c\ 
J ^ -j==p dz < max(l, |c|) 



3.4. Proofs. Fix an integer m > 2 and let Zi, 1 < i < m, be (not necessarily 
distinct) points in the complex plane that do not coalesce into a single one. Given 
i~i > 0, 1 < i < m, such that X)»=i i~i = 1 we define a meromorphic function 

rn 
i— 1 

Functions of this type are sometimes called generalized derivatives in the sense of 
Sz.-Nagy (see [3]) and may be interpreted as the resulting electrostatic force of a 
planar charge configuration (cf. loc. cit.). One of the key ingredients in our proofs 
is [3l Theorem 2], which we restate as follows: 

Lemma 1. Let <E> be as in (|3 . 7[) and denote its zeros by Wj, 1 < j < m— 1, where 
it is understood that Zi counts as a "zero" of $ of multiplicity mi — 1 if it occurs 
precisely mi times in (|3.7p . Then 



m—l 

E 

i=i 



f{wj) < - Ti)f(zi) (3.8) 



for any convex function f : C 



Remark 10. The arguments in [3] further imply that if / is a strictly convex function 
such that equality is attained in (|3.8[) then the Zj's must be collinear. 



We emphasize an important special case of Lemma [TJ 
Corollary 7. If P £ C[z] is such that degP = d > 2 i/ien 

(d-i + 1) fH<(d~i) E /(«) 

/or any convex function f : C — ► R and 1 < i < d — 1 . 

To prove Theorem 0] let (S(z), V(z)) be a pair of (higher) spectral polynomials 
of d(z) as in §3.21 set 

Z(S®) = (*«•• .,-&), l<<<n-l, 
and note that by ()1.2j) equation (|1.4p may be rewritten as 

g(g)g(f) 5W(z) 01(g) = "^ 1 1 A « ; 

Qa(z)S(*- 1 )(z) S(^)(z) + Q 2 (z) z-sf- l) 

Since Z{VS) = Z{V) V 2(5) (by the convention made at the beginning of @ and 
o-i > 0, 1 < I < p, we deduce from Lemma[T]that for any convex function / : C — > R 
one has 

r n / i \ ™— fe+l P / \ 

E«+E/(^ K E /(^n+Efi-^)^)- (3.io) 

i=l j=l v fc/ j=l i=i v fc/ 
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On the other hand, by Corollary [7] we know that 

n— fc+1 n 

n f( s t 1] ) <(n-fc + l)^/( Si ), 
i=i 3=1 
which combined with (|3.10p yields (|3.5p after some straightforward computations. 
Theorem [3] follows from the above simply by letting k = 2. 

Since functions of the type C 9 w i— ► |u>| m , m 6 Z+, and conv(Z(Q2)) 9 w 
— log|z — w\, z ^ conv(Z(Q2)), are convex, Corollaries [1] and 2] are immediate 
consequences of Theorems [3] and HI respectively. Corollaries [5] and \5\ follow from 
Corollaries [T] and HI respectively, by using Remark[3]and the expression for a k (i.e., 
afe = a(n,p, k)=n~k + l + J>2i=i a i) an( ^ noticing that 

(p — k)a>k(n,p,k) p — k (k — l)ctk(n,p, k) + n — 1 k 

and 



(p — l)ak{n,p,k) +n — 1 p (p-l)ctk(n,p,k) + n-l p 

as rt — > oo, p being fixed. To prove Corollary E] (hence also Corollary [21) we use 
again Corollary [4] and Remark [3] together with the fact that for any r C N such 

that each of the sequences < u v > , \ u„ > and < u„ > weak-* 
converges the following holds: 

,. (p- k)a k (n,p,k) 

lim r TT — = 1. 

rap^oo [p - l)a>k{n,p, k) + n(p) - 1 

k — 1 (p- k)ak(n,p, k) + n(p) — lk 

— 7 \ 1 T\ T~\ — — n P — 

p-1 {p — l)ak(n,p,k)+n{p) — l p 

M Q2 , p (/) " M Q2 . P (/) | < L p[(p-i) afc („,p >fc ) +n (p)-i] 3g 

2 

< -max 1/(^)1, 

p — 1 zeK 

where K is a compact subset of C as in Corollary El / : K — > R is any (con- 
tinuous) convex function, /i Q2 denotes the root-counting measure of Q 2 . P and 
p, Q2 p S -M(-ftT) is defined by 

supp(/i Q2p ) = Z(Q 2 , P ) = {0,p}f=l: 
ak(n,p, k) - 04 
(p- l)a k (n,p, k) +n(p) - 1 

Turning to the proof of Theorem recall first (cf., e.g., [2^]) that Ph" 1 '^ satisfies 
Jacobi's equation, i.e., the homogeneous second order linear differential equation 

(1 - z 2 )y"(z) + [(3 -a- {a + (3 + 2)z]y'(z) + n(n + a + (3 + l)y(z) = 0. 
Therefore, if Z^Pt'^') = {C n ,j}]=i denotes the zero set of P^' then 

n(n + a + (3 + l)Pt' P) (z) _ P^ ] " (z) (3 - a - (a + (3 + 2)z 



Mq 2 _ p ({0,p}) = 7Z — r„ „ i.\ I —\ 7' 1 < 1 <P- 



(l- Z 2 )Pi a ' 0y (z) Pt^'iz) l-^ 2 

1 a + 1 (3+1 



Since a > — 1 and /3 > — 1 we may apply Lemma [T] to get 

n n—1 

(n+a+(3+l)J2f({n,i) < (n+a+/3)X)/(Cj)+(»+^)/(l)+(n+a)/(-l) ( 3 - n ) 

1=1 j=l 
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for any (continuous) convex function / : [— 1, 1] — * K (cf. Remark [2] in @. Now by 
Corollary [7] we know that 

n — 1 n 

«£/(£,;,■)< fa- i)£/(Gm), 

3=1 i=l 

which together with p. lip immediately gives the desired conclusion. 

4. Further Results and Related Problems 

1. One can actually obtain stronger albeit somewhat less transparent versions of 
Theorems [3HH and establish convex domination relations for vectors whose coordi- 
nates are symmetric functions on (subsets of) Z(V)V Z(S) and Z{S t ^ k ~ 1 " 1 ) V Z(Q 2 ), 
respectively. Indeed, given deN and e S Z + with e < diet Hd, e denote the e-th 
elementary symmetric function on d elements. Then [31 Corollary 3] shows that the 
zeros and poles of the function $ defined in ()3.7p satisfy the following inequalities. 

Lemma 2. Let $ be as in (|3.7p and denote its zeros by Wj, 1 < j < m — 1, where 
as before it is understood that Zi counts as a "zero " of $ of multiplicity — 1 if it 
occurs precisely rrii times in (|3.7p . Then 



E 

l<ji<—<3d<m— 1 



/(n <e (w J1 ,...,w 3 j) 




1 ~E T " /( n ^e(2 4l ,---,^J) 



;=i 



/or any convex function f : C — > R, d € {1, . . . , m — 1} and e € Z + , e < d 



Note that Lemma Q] corresponds to the case when e = d = 1 in Lemma [21 Using 
the latter with m = n — k + p + 1 and 

{a(n,p, fc) -1 , 1 < J < — fc + 1, 
a(n,p, k)~ 1 a 3 ~ n+ k-i, n — k + 2 < j < n — k + p + 1, 

{zi\Zi +p+1 = 2(5< fe - 1 )) v z(g 2 ), {^•}^- 1 fe+p = V Z(5) 

together with (|3.9p one can then deduce weighted majorization relations of the 
aforementioned type that strengthen (|3.10p in various ways. For special choices of 
the function / one can slightly simplify these relations by first separating elements 
in Z{S^ k ^) from those in Z{Q2) (once the terms occurring in the right-hand side 
of the above inequality are appropriately regrouped) and then using Corollary [7] 
in order to estimate from above all resulting expressions that contain elements in 
Z{S^ k ~^) by means of similar expressions involving only elements in Z(S). Such 
simplifications can be made e.g. for multiplicative convex functions of the form 
f(z) = \z\", q G Z+. 

2. In view of the above results it would be interesting to know whether similar 
properties with respect to the Choquet order also hold for spectral polynomials of 
more general classes of (Lame-like) operators. Let 

*>(*) = £Q<(*)^ (4-1) 

rn 

be a linear ordinary differential operator of order k with polynomial coefficients. 
Following the terminology that we already employed for (|1.3p (cf. [TU]) we call 0(z) a 
higher Lame operator if its Fuchs index r :— max m <i<fe(deg Qi — i) is non-negative. 
If r = then d(z) is usually referred to as an exactly solvable operator in the physics 
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literature. A higher Lame operator D(z) given by (|4.1[) is said to be non- degenerate 
if deg Qk = k + r, which is equivalent to the (quite natural) requirement that 
d(z) has either a regular or a regular singular point at oo. For such an operator 
one may then consider the multiparameter spectral problem stated in (|1.4j) and 
the corresponding notions of (n-)solvability and higher Lame (i.e., Van Vleck and 
Heine-Stieltjes) polynomials. A systematic study of the latter was recently made 
in [TO]. In particular, in op. cit. it was shown that a non-degenerate higher Lame 
operator is n-solvable for all sufficiently large n and it was further proved that 
if the coefficients of D(z) are algebraically independent then for any n £ N there 
are exactly (™^ r ) Van Vleck polynomials and as many degree n Heine-Stieltjes 
polynomials, thus generalizing Heine's result (cf. Remarks [5] and [7]) . 

Problem 1. Extend Theorems [3] and H] to non-degenerate higher Lame operators 
(subject to appropriate conditions). 

An important class of linear operators which seems particularly well-suited for 
Problem [1] consists of non-degenerate higher Lame operators that also preserve 
hyperbolicity (HPOs). Indeed, as we already mentioned in Remark [7] a complete 
classification of all HPOs - i.e., linear operators T on M.[z] such that T(P(z)) has 
all real zeros whenever P € M[z] has all real zeros - was recently obtained in [8]. 
Moreover, various properties and characterizations of HPOs that belong to the 
Weyl algebra Ai (that is, operators of the form l|4.ip ) were established in [5]. For 
instance, in op. cit. it was shown that the coefficients Qi(z) of such an operator have 
all real zeros and satisfy interlacing properties like those in (|1.2|l . (Note e.g. that 
if Q2(z), Qi(z) are as in (|I.2jl and Z{Q-i) C R then the corresponding operator 
D(z) given by (| 1 .3|> is an HPO.) Furthermore, in [IU] it was proved that if a non- 
degenerate higher Lame operator d(z) is also an HPO then D(z) is solvable and 
all its Van Vleck and Heine-Stieltjes polynomials have simple real zeros. Finally, 
[H Conjecture 1] claims that HPOs either preserve or reverse the Choquet order 
on real-zero polynomials and that in particular, if T)(z) is an HPO of the form 
dUl} then Z(d{z)P(z)) -< Z(D(z)Q(z)) whenever P, Q e R[z], degP = degQ, 
Z(P), Z(Q) C R, Z(P) -< Z(Q). For results supporting this conjecture, see HOGS. 
Problem [1] should therefore be particularly interesting for non-degenerate higher 
Lame operators of HPO type. 

3. Let P n (z), n S Z+, be polynomials orthogonal with respect to a weight function 
u> supported on a (finite or infinite) interval [a, b] with u>(z) > 0, z G (a, b). It is 
well known that such polynomial families satisfy a 3-term recurrence relation 

zP n (z) = a n+1 P n+1 (z) + b n P n (z) + a„P n _i(z) 

and that under some mild assumptions on the weight function ix> (see, e.g., |19j ) 
they also satisfy a differential recurrence relation of the type 

P' n {z) - A n (z)P n ^{z) - B n {z)P n {z) 

and a second order differential equation of the form 

Plliz) + R n (z)P^{z) + S n {z)P n (z) = 0, 

where {A„(z)}„ eZ+ , {B„(z)}„ eZ+ , {i?„(z)} IlGZ+ , {S n (z)} nf zz + are certain function 
sequences. It is therefore natural to ask the following questions. 

Problem 2. Investigate whether there are weighted majorization relations similar 
to those in Theorem [S] and /or weighted majorization relations involving the zeros 
of any two (or three) consecutive terms for 

(i) Laguerre-like and Hermite-like polynomials; 

(ii) (appropriate classes of) general orthogonal polynomials. 
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